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Abstract 

We construct extensions of the Standard Model with two Higgs 
doublets, where there are flavour changing neutral currents both in the 
quark and leptonic sectors, with their strength fixed by the fermion 
mixing matrices Vckm and Vpmns- These models are an extension 
to the leptonic sector of the class of models previously considered by 
Branco, Grimus and Lavoura, for the quark sector. We consider both 
the cases of Dirac and Major ana neutrinos and identify the minimal 
discrete symmetry required in order to implement the models in a 
natural way. 
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1 Introduction 



Understanding the mechanism of gauge symmetry breaking is one of the 
fundamental open questions in Particle Physics. Indeed, even if elementary 
scalar doublets are responsible for the breaking of the gauge symmetry, one 
does not know whether the breaking is generated by one, two or more scalar 
doublets. This question is specially relevant in the beginning of the LHC 
era, with the prospects of experimentally probing the mechanism of gauge 
symmetry breaking. 

In the Standard Model (SM) there is only one Higgs doublet and, as a 
result, scalar couplings are automatically flavour diagonal in the quark mass 
eigenstate basis. In multi-Higgs models, this is no longer true and one is 
confronted with dangerous flavour changing neutral currents (FCNC) at tree 
level, unless one introduces a symmetry leading to natural flavour conserva- 
tion (NFC) [1], [2], [3], in the Higgs sector or some alternative mechanism to 
naturally suppress FCNC. 

A possible alternative scenario for suppressing FCNC is through the as- 
sumption that all flavour violating neutral couplings be proportional [5] 
to small entries of the Cabibbo-Kobayashi-Maskawa matrix (Vckm) [6], [7]. 
This is one of the ingredients of the Minimal Flavour Violation (MFV) prin- 
ciple, introduced for the quark sector in Refs. [S], [H] and later on extended 
to the leptonic sector [TU], [TT], [12]. The first models of the MFV type, 
in the framework of two-Higgs doublets and without ad hoc assumptions, 
were proposed by Branco, Grimus and Lavoura (BGL) [13]. In the BGL 
models, the MFV character results from an exact discrete symmetry of the 
Lagrangian, spontaneously broken by the vacuum. Another proposal for the 
structure of the scalar couplings to fermions is the suggestion that the two 
Yukawa matrices are aligned in flavour space |14] . 

Recently, we have proposed an extension of the hypothesis of MFV to 
general multi-Higgs models with special emphasis on two Higgs doublets [T5] . 
In that work there is a detailed analysis of the conditions for the neutral 
Higgs couplings to be only functions of Vckm elements as well as a MFV 
expansion for the neutral Higgs couplings to fermions. This expansion is 
built by combining the most basic elements [16j that transform appropriately 
under weak basis transformations, with terms proportional to the fermion 
mass matrices. 

In this paper, we study how our analysis can be extended to the leptonic 
sector, considering both the case where neutrinos are Dirac particles and the 
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case where neutrinos are Majorana particles, acquiring naturally small masses 
through the seesaw mechanism. Note that this extension to the leptonic 
sector is crucial in order to study stability under renormalization as well as 
to do a full phenomenological analysis of BGL type models. 

This paper is organized as follows: In section 2 we extend the BGL model 
to the leptonic sector with the imposition of lepton number conservation. 
Therefore, in this case, neutrinos are Dirac particles. Furthermore, we show 
that in this case the one loop renormalization equations for the corresponding 
Yukawa couplings, both in the quark and leptonic sector, obey the equations 
that guarantee the dependence of Higgs FCNC solely on functions of the 
mixing matices. In this section we also deal with the question of the unique- 
ness of BGL models. In section 3 we extend the BGL model to the leptonic 
sector without the imposition of lepton number conservation. We start by 
discussing the effective low energy scenario with Majorana neutrinos, and its 
stability. Next, we analyse the leptonic sector in the seesaw framework taking 
into consideration both the low and high energy couplings to the Higgs fields 
involving neutrinos. In section 4 we argue in favour of having the symmetry 
leading to MFV softly broken in the scalar potential. Finally, in section 5 
we present our Conclusions. 

2 Minimal Flavour Violation with Dirac Neu- 
trinos 

2.1 Framework 

The extension of BGL models to the leptonic sector depends on the neutrino 
character. In this section, we analyse the leptonic sector of models that 
account for neutrino masses by enlarging the Standard Model (SM) through 
the introduction of three righthanded neutrinos while at the same time 
imposing total lepton number conservation. As a result, only Dirac mass 
terms are generated and neutrinos are Dirac particles. We consider models 
with two Higgs doublets such that the flavour changing neutral currents are 
controlled by the Vckm matrix in the quark sector and by the Pontecorvo- 
Maki-Nakagawa-Sakata (PMNS) matrix [T7], [IB], [IH] in the leptonic sector. 
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The full Yukawa couplings are given by: 

ni$i/0 - U[ U2^2l% - ^ Si$iZ/0 - S2$24 + h.C. ,(1) 

where Fj, Aj denote the Yukawa couplings of the lefthanded quark doublets 
Ql to the righthanded quarks (i^j, and to the Higgs doublets Ilj, Ej de- 
note the couplings of the lefthanded leptonic doublets L° to the righthanded 
charged leptons Z^, neutrinos and to the Higgs doublets. Lepton num- 
ber conservation prevents the existence of invariant mass terms of Majorana 
type for righthanded neutrinos. These will appear in section 3 in the see- 
saw framework, where the requirement of lepton number conservation will 
be relaxed. 

In order to obtain a structure for Fj, Aj such that there are FCNC at 
tree level with strength completely controlled by Vckm, Branco, Grimus and 
Lavoura imposed the following symmetry on the quark and scalar sector of 
the Lagrangian [IB] : 

Qlj exp (ia) Q\ , u\ exp {i1a)u\ , $2 ^ exp (ia)$2 , (2) 

where a 7^ 0, tt, with all other quark fields transforming trivially under the 
symmetry. The index j can be fixed as either 1, 2 or 3. Alternatively the 
symmetry may be chosen as: 

exp {ia) Qlj , d^j^j exp {i2a)d%^^ , $2 exp {-ia)^2 ■ 

(3) 

The symmetry given by Eq. (|2]) leads to Higgs FCNC in the down sector, 
whereas the symmetry specified by Eq. ^ leads to Higgs FCNC in the up 
sector. The neutral Higgs interactions with the fermions, obtained from the 
quark sector of Eq. ([1]) are given by 

£y (neutral, quark) = -IFJ- [MdH^ + N^R + iN^J] 4 + 

- <^[M„/7° + iV°i? + *iV°/]4 + h.c. , (4) 

where v = \/vf + V2 = {V^Gf)~^^'^ ~ 246 GeV, Gp is the Fermi coupling 
constant and if °, R are orthogonal combinations of the fields pj, arising when 
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one expands [20] the neutral scalar fields around their vacuum expectation 
values (vevs), 0° = ^{vj + pj + ir]j), choosing in such a way that it 
has couplings to the quarks which are proportional to the mass matrices, as 
can be seen from Eq. (jl]). Similarly, / denotes the linear combination of rjj 
orthogonal to the neutral Goldstone boson. The mass matrices and 
and the matrices and A^^ are given by: 

M, = -^{viTi + v^e^'T^) , M„ = ^^Ai + v^e-^'A^) , (5) 

here 6 denotes the relative phase of the vevs of the neutral components of 
The matrices M^, are diagonalized by the usual bi-unitary transfor- 
mations: 

Ul^MdUdR = diag (m^, m^, nib) , (7) 

Ul^MJJuR = Du = diag (m„, mc, rrit) . (8) 

The flavour changing neutral currents are controlled by the matrices Nd and 
Nu related to A*"^ and by the following transformations: 

Nd = Ul^N'dUm , iV„ = Ul^KU^R . (9) 

In the case of the symmetry given by Eq. ([2]), for j = 3 there are FCNC in 
the down sector controlled by the matrix A^^ given by [13] 



(A^,)., ^ ^{DdU - ( ^ + ^ ) iVij,j,MVcKM)3ADdh ■ (10) 

whereas, there are no FCNC in the up sector and the coupling matrix of the 
up quarks to the R and / fields is of the form: 

Nu = -—diag (0, 0, nit) + —diag (m„, ni^, 0) . (11) 

In this example, the Higgs mediated FCNC are suppressed by the third row 
of the VcKM matrix, therefore obeying to the additional constraint imposed 
on models designated as of the MFV type [21]. 
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As shown in reference [12], the symmetries given by Eq. ([2]) or by Eq. 
lead to 

v]r, = r,, v]r, = o, (12) 

r]A, = A2, P]A, = 0, (13) 

where 7 stands for u (up) or d (down) quarks, and Vj are the projection 
operators defined by 

VJ = U^lP.uI, , Vf = U,,P,Ul , (14) 

and {Pj)ij^ = Sjidjk- Note that Eqs ffT^ and f[T51) . guarantee that the Higgs 
fiavour changing neutral couplings can be written in terms of quark masses 
and VcKM entries [Ldj. This is a crucial feature for BGL models to be con- 
sidered as of Minimal Flavour Violation type (MFV). 

In the leptonic sector, with Dirac type neutrinos, there is perfect analogy 
with the quark sector, consequently MFV is enforced by one of the following 
symmetries. Either 

Lik exp (ia) L% , z/^;. exp {i2a)h'%^, , $2 ^ exp (za)<l>2 , (15) 

or 

Llk exp {ia) L% , 1% exp {i2a)l% , $2 ^ exp (-ia)$2 , 

(16) 

where, once again, a 7^ 0, vr, with all other leptonic fields transforming triv- 
ially under the symmetry. The index k can be fixed as either 1, 2 or 3. 
Similarly, for the leptonic sector, these symmetries imply 

pf n2 = 1I2 , pf Hi = , (17) 

PfS2 = S2, PfSi = 0, (18) 

where P stands for charged lepton or neutrino. In this case 

Vi = UiLPkUl , VI = U.LPkUU , (19) 

where U^l and Uil are the unitary matrices that diagonalize the correspond- 
ing square mass matrices 

Ul^MiMjUiL = diag {mlmlml) , 
Ul^M^MlU^L = diag{ml^,ml„ml,) , (20) 
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with Ml and M„ of the form 



2.2 Renormalization Group Study 



(21) 



Equations ( !T2|) and (fT3|) together with Eqs. ( !T7|) and ( !T8|l guarantee that the 
Higgs FCNC are functions of fermion masses and of the CKM and PMNS 
matrices. Therefore, it is crucial to guarantee the stabihty of these equations 
under renormahzation. 

The one loop renormalization group equations (RGE) for our Yukawa 
couplings can be generalized from reference [22] to 



arTk + 

2 



+ J2 [sTr (r,rt + At A,) + Tr (n^n} + St 



1=1 

2 



r/ + 



+ (-2A, At r, + r.rtr, + ^ a, Atr, + ^r.rjr, j 



(22) 



VAk = oaA^ + 



+ J2 pr (a, At + rt r,) + Tr (s,st + nt n, 
1=1 

+ J2 (-2rirt Ai + Afe AtA, + ^r^r} A, + ^a^aJa^^ , (23) 
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+ [3Tr(r,rt + AtA,) +Tr(n,nt + sts, 

1=1 

+ (-2S'^Inz + n^ntn^ + ^s^stn^ + ^n^ntn^") , (24) 
1=1 ^ ^ 
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+ J2 [sTr ( A, Aj + rir,) + Tr (s^sf + nln, 

1=1 



where 15 = IGvr^/i (d/dfi) and /i is the renormahzation scale. The coefficients 
ar, ca, On and as are given by [23] : 

an = -^9 , os = -^9 - -^9 , (26) 

where 5^^, 9 and g' are the gauge couphng constants of SU{3)c, SU{2)l and 
U{1)y respectively. To show that Eqs. ([H]), ([13]), ([17]) and ([H]) are stable 
under RGE one has to show that 

V] {VT2) = {VT2) , V] (OTi) = , (27) 
P;(PA2) = (I5A2) , P7(PAi) = 0, (28) 

Pf (PHs) = {VIi2) , Pf (PHi) = , (29) 
Pf (DE2) = (PS2) , Pf (PEi) = , (30) 

which guarantee that the Yukawa couplings at each different scale still verify 
equations of the same form. 

It is interesting to notice that if one does not use the conditions for the 
leptonic sector given by Eqs. f[T7]) and f lTS]) one is lead, for example, to: 

V] (PEi) = Tr {tI^U\ + Sl S2) , (31) 

also, one must use the equality: 

Trfni^ + SjSa) =0 (32) 



in order to show that P" (Pr2) = {^2)- Clearly, Eq. fl32]) follows from 
Eqs. f lTTj) and f[TH]) . since in this case we have: 



Tr(^nin5 + sls2) = Tr(^nin^p/ + slpfs2^ 

= Tr (P/Hi) + (P/Si) ^ S2^ = , (33) 



so that Eq. f l32p is enforced by the MFV leptonic conditions. It is the entire 
set of equations both in the quark and in the leptonic sector that guarantee 
the stabihty of these models. This fact should not come as a surprise since the 
relations given by Eqs. f|T2|) . f lT3|) . f|T7|) and f|T8|) follow from the imposition 
of a symmetry on the full Lagrangian. 

In the quark sector there were six possible different implementations of 
BGL type models. Three with FCNC in the down sector, each one corre- 
sponding to a different choice for the index j, and three with FCNC in the 
up sector also for the different choices of j. In the extension to the leptonic 
sector with Dirac neutrinos one has another set of six different leptonic im- 
plementations obtained in a similar fashion. In total, one may consider thirty 
six different MFV models of BGL type in the case of Dirac neutrinos. 

In Ref . |15] we presented a MFV expansion for and built with terms 
proportional to Ma and respectively, as well as products of terms which 
transform like H^i and under weak basis transformations, multiplying 
and Mu- We identified Vj and as the simplest such elements with the 
appropriate transformation under changes of weak basis. In fact, and 
can be decomposed as [16j : 

Hdiu) = J2<iu),Pf^''^ ■ (34) 

i 

As a result we obtained, for example, simple models of MFV type with 
Higgs mediated FCNC in both sectors, like the one given by the following 
equations: 

iV,° = ^ Md- (- + -] UuLmi^ Md , (35) 
iV° = ^ M„ - (^^ + UdLPiUl^ Mu . (36) 

Several different possible variations beyond BGL models were considered in 
Ref. [U], obtained from different combinations of and with and 
Mu- We pointed out in Ref. [l^ that the zero texture structure of these 
models is more involved than in the BGL case and that the question of as- 
suring its loop stability, through the introduction of symmetries, was not 
obvious. We can now address this question with the help of the renormal- 
ization group equations for the Yukawa couplings given by Eqs. 
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and fl25l) . All these additional MFV models, as well as the BGL models, 
lead to relations of the following form: 

Is = T2 , Fi = , (37) 

PfA2 = A2, PfAi = 0, (38) 

and similar equations for the leptonic sector. For a = /3 and i = j we are in 
a BGL model in the quark sector. Models with a 7^ /3 or i 7^ j correspond 
to additional cases presented in Ref. [T3] . 

It can be readily verified that, in general 

(OTi) = -^PfAiAlFi-2PfpfA2AlPfF2 + ^PrpfA2A5pfFi . (39) 

We have already shown that for BGL models we have 

(PFi) = . (40) 

In the case a = f3, i j we have 

vrv^ = , (41) 

due to the fact that these are projection operators. So we are left with 

Vr (OTi) = -^PfAiAlFi , (42) 

which, in general, is different from zero. Therefore we conclude that this type 
of models cannot be enforced by symmetries. The consideration of equation 
Vf {T>T2) = ('CF2) would lead to similar difficulties and therefore, would 
allow us to draw a similar conclusion. As a result we may conclude that out 
of the models described by Eqs. fl37|l and f p8|) and their generalization to the 
leptonic sector, only BGL type models can be enforced by some symmetry. 
The same question was recently addressed in Ref. [24J following a different 
approach. There it was shown that BGL models are the only ones that 
survive among a large set of models enforced by abelian symmetries. 
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3 Minimal Flavour Violation with Majorana 
Neutrinos 



3.1 Low Energy Effective Theory and Stability 

In the previous section, we assume that neutrinos are Dirac particles. An 
alternative possibility is to allow for lepton nonconservation leading to an 
effective Majorana mass term for the three light neutrinos of the form 

1 T 

'^Majorana = C'^rn^ul + h.c. , (43) 

which violates lepton number. Such a mass term is generated after sponta- 
neous gauge symmetry breaking from an effective dimension five operator O 
which, in the two Higgs doublet model can be written as: 

0=Y. E E {LLa4'^C~'L^Pc){e'^%,){e^%,). (44) 

ij=l a,f3=e,fi,T a,b,c,d=l 

This operator contains two lefthanded lepton doublets and two Higgs dou- 
blets and can be viewed, for example, as arising from the seesaw mechanism 
after integrating out the heavy degrees of freedom. In the seesaw context 
the heavy degrees of freedom are the righthanded neutrinos. The seesaw 
framework will be analysed in the next subsection. 

In this context we have, in the leptonic sector, the two flavour structures 
introduced before: 

Cy, = -iJi ni<i>i/^ - U[ n2$2/^ + h.c. , (45) 

together with the four new flavour structures given by the k*^*-^-' matrices. 
A priori, it looks more difficult to implement MFV in the case of Majorana 
neutrinos. However, this can be done by imposing the following Z4 symmetry 
in the effective Lagrangian including the terms given by Eqs. (jH]) and fH5|) : 

Llj -)> exp (ia) , $2 -> exp (m)<l>2 , (46) 

with a = 71 /2. Imposing this Z4 symmetry implies: 









(47) 
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and taking for definiteness j = 3 we get 



K 



(11) 



X X 
X X 




K 



(22) 





X 



(48) 



(22) 

fixing the angle a as 7r/2 ensures that 7^ so that the determinant 
of the resulting neutrino mass matrix does not vanish automatically. The 
Majorana mass matrix for the neutrinos is given by: 



1 
2 



2 ^ 2 ^ 



This ^4 symmetry also implies the following structure for IIi and 112: 



Hi 



X 
X 





X 
X 





X 
X 





Ho 






X 






X 






X 



(49) 



(50) 



The neutrino mass matrix rrij, is block diagonal with each block given by 
a different n matrix. As a consequence, in the diagonalization of m^, the 
matrices k^^^^ and k*^^^^ are diagonalized separately. Therefore, any linear 
combination of these two matrices will be simultaneously diagonalized. As a 
result the lepton number violating Weinberg operator [25] of Eq. f H^ does 
not give rise to Higgs mediated FCNC in the neutrino sector. For the charged 
lepton sector the situation is similar to the one encountered in the previous 
section for the symmetry given by Eq. ( |T5|) . leading to Higgs mediated FCNC 
in this sector. 

The symmetry imposed by Eq. ( H6|) in the effective low energy theory 
leads, for j = 3 for instance, to the following conditions: 



,(12) 



K 



(21) 











V!(ni = 



/^(22)p^^' 



K 



(22) 



(51) 



It can be easily verified, as we have done in section 2.2, and following the 
RGE presented in Ref. [26] that these equations are indeed stable under 
renormalization, since they keep the same form at all scales. 



3.2 Seesaw Framework 

In this section, we analyse the leptonic sector in the seesaw framework [27], 
[25] , [2n] , [30] , [SI] ■ We include one righthanded neutrino per generation and 
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do not impose lepton number conservation. The leptonic part of Yukawa 
couplings and invariant mass terms can then be written: 

+ l^^'/c-'Mniy'R + h.c. . (52) 

The matrix Mr stands for the righthanded neutrino Majorana mass matrix. 
The leptonic mass matrices generated after spontaneous gauge symmetry 
breaking are given by: 

mi = -^{viUi + V2e''U2) , nio = ^(^i^i + V2e-''^2) • (53) 

Note that the notation has changed from the one in section 2, we now have 
mi = Ml and m £> replaces Mj, in order to avoid confusion with light neutrino 
masses in the seesaw framework. The leptonic mass terms obtained from 
Eq. (!52|) can be written as: 



£ma8s = -ll m, 1% + , (Af) C-'Af ' ( ) + h.c. . (54) 



M = I °r 1 . (55) 



with 

, mi Mr 
We use the following convention: 

{^lY ^ C^^i^iT ■ (56) 
The charged current couplings are given by: 

Cw = -^W;lf^r 1^1 + ^.0.. (57) 

The neutral Higgs interactions with the fermions, obtained from Eq. flS2]) can 
be written: 

£y (neutral, lepton) = -/« ^ [miH° + Nj>R + iN^I] /° + 

- [moH^ + KR + iN^J] vl + h.c. , (58) 
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with 

iV° = ^n, - ^e^^n, , (59) 
iV° = ^S, - ^e-^S, . (60) 

There is a new feature in the seesaw framework due to the fact that in the 
neutrino sector the hght neutrino masses are not obtained from the diago- 
nahzation of ttid- 

The 6x6 neutrino mass matrix Ai is diagonahzed by the transformation: 

V^M*V = 2i , (61) 

where 3) = dia.g{m^i,m^2^^u3, Mi, M2, M^), with m^i and Mj denoting the 
masses of the physical hght and heavy Majorana neutrinos, respectively. It 
is convenient to write the matrices V and ^ in the following block form: 



^=(s rj • ®=(o ■ 

In the seesaw framework, with the scale of ^ v the matrix K coincides 
to an excellent approximation with the unitary matrix U that diagonalizes 
the effective mass matrix rrieff for the light neutrinos: 

rrieff U* = d with f^eff = — "^dyit^^ • (^'^) 



The matrix G verifies the exact relation 

G = mDT*D-^ , (64) 

while S is given to an excellent approximation by [32] : 

= -K^moM^^ , (65) 

It is clear from Eqs. (Ell) and (l65ll that G and 5* are of order rriD/MR, therefore 
strongly suppressed. This in turn means that the 3x3 matrices K and 
T are unitary to an excellent approximation. The matrix T is also very 
approximately determined by: 

T^MrT* = D . (66) 
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The physical fermion fields /, v and are then related to the weak basis 
fields by: 



— ViiIl , 1% — Uijil 



Uul + GNl, vl = S*vl + T*Nl. (67) 



In terms of physical fields the charged gauge current interactions become 



w 



9 ,T- 



-^{Il^^U^ulW^^ + IlI.QNlW) + h.c. . 



(68) 



Uu denotes the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix defined 
by the product (Ui^^U). The second term of Cw in Eq- (EH]), with mixing 
given by Q = (t/^I^G) is suppressed by G and involves the heavy neutrinos 
N which are not relevant for low energy physics. 

In general the couplings of Eq. f l58|) lead to arbitrary scalar FCNC at tree 
level. In order for these couplings to be completely controlled by the PMNS 
matrix we introduce the following Z4 symmetry on the Lagrangian: 

L% exp (ia) L^^ , ^ exp {i2a)v%^^ , $2 ^ exp {ia)^2 , (69) 

with a = 7r/2 and all other fields transforming trivially under Z4. The most 
general matrices Ilj, Sj and Mr consistent with this Z4 symmetry have the 
following structure: 





X 


X 


X 




■ 





Hi = 


X 


X 


X 


n2 = 





















X 


X 




X 


X 


" 




■ 





Si = 


X 


X 





S2 = 
























X 








(70) 



R 



X X 
X X 
X 



1) 



where x denotes an arbitrary entry while the zeros are imposed by the sym- 
metry Z4. Note that the choice of Z4 is crucial in order to guarantee M33 7^ 
and thus a non- vanishing det Mr. The same choice was required in the pre- 
vious subsection in order to allow for a non-vanishing determinant for the 
effective Majorana neutrino mass matrix. In this weak basis the following 
important relations are verified: 

" 

PaHs = , PaHi = , with P3 









as well as 



AS 



3^2 



P3S1 = 



(72) 



(73) 
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3.3 Full Seesaw Higgs couplings 

Let us now write the neutral scalar couplings of the charged leptons in the 
mass eigenstate basis: 

H° - 

£y (neutral) = / Di I 

- ^ I {NnR + nI^l)) l + t^-l {Nhr - iV/7L)) / , (74) 

where 7l = (1 - 75)/2 ,lR=il + 75)/2 and Ni = Ui{ iV" Uir. 

The fact that U given by Eq. ( l63i) is block diagonal with no mixing for 
the third family leads to: 

(iV,),, ^ {U,i U,nh = -(Dih '(- + -) {UlUU^hiDih ■ (75) 

Uy is the PMNS matrix. 

We obtain the neutrino couplings to neutral scalars from the last term of 
Eq. (!58|) . It is useful to rewrite A^^ in the form: 

K = —^D - ^ — + — e S2 . 76 

Notice that the first term of A^^ is proportional to rrio and therefore these 
couplings to the fields R and / have a structure similar to the couplings 
in Eq. (1551) . The couplings of the neutrino mass eigenstates z/j (light), 
(heavy) to the neutral scalars H^, R and I are more involved than the cou- 
plings of the charged leptons, since they include light-light, light-heavy and 
heavy-heavy couplings. In the sequel, we shall consider each one of these 
terms, displaying their explicit form in the present model. 

couplings 

(i) light-light couplings. 
These couplings can be written 

= ^ ViZH'^^, + h.c. , (77) 

V 

where 

A = U^moS* = d . (78) 
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These couplings among light neutrinos are flavour diagonal and are propor- 
tional to the light neutrino masses. From the point of view of the effective 
low energy theory there are no scalar FCNC in the neutrino sector, since, 
as will be shown, the term of in ZI2 given by Eq. f l76|) . corresponding to 
light-light couplings {WT12S*) will not generate nonzero off-diagonal entries. 
Its form is given exphcitly, in the sequel, by Eq. fl57|) . 



(ii) light-heavy couplings. 
We write these terms as 

B 



J y J 

E = G^moS* 



where 



From Eqs. 



h.c. 



B = U^moT* , , 
and fl66p one can write 

B = (iVdO 



D) 



(79) 



(80) 



(81) 



where is an orthogonal complex matrix. This expression readily follows 
from the Casas and Ibarra parametrization [33]. The fact that as well 
as Mr are block diagonal implies that is also block diagonal and can be 
parametrized as: 

cos Z ± sin Z 
0^= -sinZ ±cosZ , (82) 
1 

with Z complex. These couplings, given by the matrix B, are not suppressed 
by the mixing matrices but the fact that the heavy neutrino fields, N , have 
masses of order Mr implies that they cannot be produced at low energies. 
Using Eqs. ( 16^ and ( l65l) it can be readily verified that: 



E = -D-^iVdO'VoYd 



(83) 



These couplings, given by the matrix E, are suppressed by both matrices G 
and S therefore they are much smaller than those given by -B, in addition, 
they also include a heavy neutrino. 



(iii) heavy-heavy couplings. 
One has for these couplings: 



.HO 



^ NT^H'Nl^ + h.c. 



M) 
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where 



F = G^ttidT* , 



(85) 



it can be readily verified that: 



(86) 



These are couphngs among heavy neutrinos and furthermore are suppressed 
by the mixing matrix G. 

R and / couplings 



Concerning the neutral couplings to R and I the first term of given by 
Eq. fl76|) leads to currents with the same structure as those mediated by 
H^. The second term of A'"^ leads to diagonal coupling matrices, due to the 
block structure of S2 given by Eq. f lTTl) and the fact that, as a result of the 
patterns given by Eq. f lTT]) for the neutrino mass matrices, the matrices U, 
G, S and T are block diagonal with no mixing in the third row and column. 
The additional couplings to R and / are derived by replacing m£) by S2 in A, 
B, E and F introduced by Eqs. ([THD, ([HOD and (fM])- From Eq. ([7HD and the 
definition of m^i given by Eq. (l53l) we obtain the following additional term 
for light-light couplings to R and /: 



From Eqs. ( IHOl) . (IHTl) and (1831) we obtain the following light-heavy coupling 



Finally from Eqs. ( l85l) . ( l86l) we obtain the following heavy-heavy coupling 




(87) 



terms: 




(88) 



(89) 



term: 




(90) 
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Combining the two contributions from A*"^, the hght-hght couphngs to R 
and / are diagonal. For the first two generations the coefficients are given by 
— — and — — respectively. For the third generation it is given by — . 

The light-heavy couplings to R and / are block diagonal. Compared to 
the corresponding couplings the block (12) is multiplied by the ratio of 
vevs ^ and the (33) coupling is multiplied by — ^. Likewise for the heavy- 
heavy couplings to R and I. 

From the point of view of low energy physics, the example given is a model 
of BGL type, with no Higgs mediated FCNC in the up sector (light neutri- 
nos) and with the strength of the FCNC in the down sector controlled by 
the PMNS mixing matrix. All flavour changing neutral couplings with heavy 
neutrinos are parametrized by both light and heavy neutrino masses, and 
the product of matrices i\fdO'^\fD, with of the form given by Eq. fl82|) . 
Heavy neutrino decays may be the source of the baryon asymmetry of the 
universe through leptogenesis [31] with sphaleron processes [35], [36] . 

Next we write the Yukawa couplings to the charged Higgs, using 
Eqs. (187]) and (|90|). 

couplings 



The charged Higgs interactions with the fermions, obtained from Eq. ( I52ll 
are given by 



£y (charged) 



{vlN^ll + vlNlhl) + h.c. . 



(91) 



V 



In the fermion mass eigenstate basis these interactions become: 



Cy (charged) 



V2H+ 



[WUlNilR + NLQ^Niln] + 



V 




(92) 
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4 The scalar Potential 



The Z4 symmetry which we have imposed on the Lagrangian, forbids vari- 
ous gauge invariant terms in the scalar potential, such as 0^02 , 0i020l0i,, 
0^020102- As a result, the Higgs potential has an exact ungauged accidental 
continuous symmetry, which is not a symmetry of the full Lagrangian. After 
spontaneous gauge symmetry breaking, the accidental symmetry leads to a 
pseudo-Goldstone boson. The simplest way of avoiding the pseudo-Goldstone 
boson is by breaking the Z4 symmetry softly through the introduction of the 
term mi24>\<p2 + h.c. This term avoids the pseudo-Goldstone boson which 
acquires a squared mass proportional to |mi2|. In order to discuss CP vi- 
olation in this class of models, one has to consider separately the cases of 
explicit and spontaneous CP violation. 

Explicit CP violation - If one does not impose CP invariance at the 
Lagrangian level, Yukawa couplings are complex. In spite of the special form 
of these couplings, due to the presence of the Z4 symmetry, it can be readily 
checked that there is in general CP violation through the Kobayashi-Maskawa 
(KM) mechanism. The simplest way of verifying that this is the case, is by 
noting that Ha = MdMl is a generic complex Hermitian matrix while Hu is 
a block diagonal matrix. One can easily compute Tr[ifu, Rdf" and show that 
in general this weak-basis invariant does not vanish thus proving [37] that 
there is CP violation through the KM mechanism. In order to check whether 
there are in this model other sources of CP violation, one has to look at the 
scalar potential. It can be readily checked that the scalar potential, by it- 
self, is CP invariant since the phase of mi2 can be removed by rephasing the 
scalar doublets, thus rendering the potential real. The powerful Higgs-basis 
invariant CP-odd conditions derived in Ref. |38] would obviously provide the 
same answer, however this is a straightforward case. In this variant of the 
model, one has all CP violation arising from the KM mechanism. However, 
note that there are, for example, new contributions to Bd - Bd apart from 
the usual box diagrams of the Standard Model. These new contributions are 
mediated by tree level scalar interactions, which are proportional to (V^feV^^)^, 
therefore with the same phase as the SM box contribution. 

Spontaneous CP violation - It can be readily checked that even in the 
presence of the soft breaking term mi2(/>l02, one cannot achieve spontaneous 
CP violation, without enlarging the scalar sector. On the other hand, one 
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may obtain spontaneous CP violation by introducing scalar singlets. However 
in order for the phase arising from the vacuum to be able to generate a 
complex CKM matrix, one has to introduce vector-like quarks [ 39] . 

5 Conclusions 

We have analysed how to extend to the leptonic sector, BGL models sat- 
isfying the minimal flavour violation (MFV) hypothesis. Both the cases of 
Dirac and Majorana neutrinos were considered. In the case of Dirac neu- 
trinos the extension to the leptonic is straightforward with great similarity 
to the quark sector. We have shown that if type-I seesaw mechanism is 
adopted, the requirement of having a non-singular Majorana mass matrix 
for the righthanded neutrinos further restricts the choice of the discrete sym- 
metry which allows for realistic BGL models in the leptonic sector. A striking 
result of our analysis is the fact that this restricted form of the symmetry is 
also required when considering the low energy effective theory with Majorana 
neutrinos. In particular, it was pointed out that BGL models satisfying the 
MFV paradigm can be extended in a natural and elegant way to the leptonic 
sector with Majorana neutrinos, through the introduction of a symmetry, 
imposed on the full Lagrangian. Furthermore we derive the equations which 
guarantee calculability of Higgs FCNC in terms of masses, Vqkm and Vpmns 
matrices showing that these equations are stable under renormalization. We 
have also analysed the scalar potential which acquires an exact ungauged 
accidental continuous symmetry arising from the absence of various terms 
forbidden by the symmetry. We have pointed out that the simplest way 
of avoiding the resulting pseudo-Goldstone boson is through the addition of a 
quadratic term in the scalar potential, thus softly breaking the Z4 symmetry. 
Finally, we emphasize that the relevance of BGL models stems in good part 
from the fact that the most general tree-level flavour violating neutral cur- 
rents are naturally suppressed by small Vckm elements like the combination 

A full analysis of BGL models is beyond the scope of this paper and will 
be presented elsewhere [10] • It is clear that the extension of BGL models to 
the leptonic sector is essential in order to make possible the above analysis 
and furthermore, to allow for a consistent analysis of the renormalization 
group evolution. 



20 



Acknowledgements 

This work was partially supported by Fundagao para a Ciencia e a Tecnolo- 
gia (FCT, Portugal) through the projects CERN/FP/109305/2009, 
PTDC/FIS/098188/2008 and CFTP-FCT Unit 777 which are partially funded 
through POCTI (FEDER), by Marie Curie Initial Training Network "UNILHC" 
PITN-GA-2009-237920, by Accion Complementaria Luso-Espanhola PORT2008- 
03 and FPA-2008-04002-E/PORTU, by European FEDER, Spanish MICINN 
under grant FPA2008-02878 and GVPROMETEO 2010-056. GCB and MNR 
are very grateful for the hospitality of Universitat de Valencia during their 
visits. FJB and MN are very grateful for the hospitality of CFTP/IST Lisbon 
during their visits. MN thanks MICINN for a Juan de la Cierva contract. 

References 

[1] S. Weinberg, Phys. Rev. Lett. 37 (1976) 657. 

[2] S. L. Glashow and S. Weinberg, Phys. Rev. D 15 (1977) 1958. 

[3] E. A. Paschos, Phys. Rev. D 15 (1977) 1966. 

[4] A. Antaramian, L. J. Hall and A. Rasin, Phys. Rev. Lett. 69 (1992) 
1871 [arXiv:hep-ph/9206205]; L. J. Hall and S. Weinberg, Phys. Rev. 
D 48 (1993) 979 [arXiv:hep-ph/9303241]. 

[5] A. S. Joshipura and S. D. Rindani, Phys. Lett. B 260 (1991) 149; 

[6] N. Cabibbo, Phys. Rev. Lett. 10 (1963) 531. 

[7] M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49 (1973) 652. 

[8] A. J. Buras, P. Gambino, M. Gorbahn, S. Jager and L. Silvestrini, 
Phys. Lett. B 500 (2001) 161 [arXiv:hcp-ph/0007085]. 

[9] G. D'Ambrosio, G. F. Giudice, G. Isidori and A. Strumia, Nucl. Phys. 
B 645 (2002) 155 [arXiv:hep-ph/0207036]. 

[10] V. Cirigliano, B. Grinstein, G. Isidori and M. B. Wise, Nucl. Phys. B 
728 (2005) 121 [arXiv:hep-ph/0507001]. 



21 



[11] S. Davidson and F. Palorini, Phys. Lett. B 642 (2006) 72 [arXiv:hep- 
ph/0607329]. 

[12] G. C. Branco, A. J. Buras, S. Jager, S. Uhlig and A. Weiler, JHEP 
0709 (2007) 004 [arXiv:hep-ph/0609067]. 

[13] G. C. Branco, W. Grimus and L. Lavoura, Phys. Lett. B 380 (1996) 
119 [arXiv:hep-ph/9601383]. 

[14] P. Tuzon and A. Pich, Acta Phys. Polon. Supp. 3 (2010) 215 
[arXiv: 1001. 0293 [hcp-ph]]. M. Jung, A. Pich and P. Tuzon, JHEP 1011 
(2010) 003 [arXiv: 1006.0470 [hep-ph]]. M. Jung, A. Pich and P. Tuzon, 
arXiv:1011.5154 [hep-ph]. 

[15] F. J. Botella, G. C. Branco and M. N. Rebelo, Phys. Lett. B 687 (2010) 
194 [arXiv:0911.1753 [hep-ph]]. 

[16] F. J. Botella, M. Nebot and O. Vives, JHEP 0601 (2006) 106 
[arXiv:hep-ph/0407349]. 

[17] B. Pontecorvo, Sov. Phys. JETP 7 (1958) 172 [Zh. Eksp. Teor. Fiz. 34 
(1957) 247]. 

[18] Z. Maki, M. Nakagawa and S. Sakata, Prog. Theor. Phys. 28 (1962) 
870. 

[19] B. Pontecorvo, Sov. Phys. JETP 26 (1968) 984 [Zh. Eksp. Teor. Fiz. 
53 (1967) 1717]. 

[20] T. D. Lee, Phys. Rev. D 8 (1973) 1226. 

[21] A. J. Buras, M. V. Carlucci, S. Gori and G. Isidori, JHEP 1010 (2010) 
009 [arXiv:1005.5310 [hep-ph]]. 

[22] P. M. Ferreira, L. Lavoura and J. P. Silva, Phys. Lett. B 688 (2010) 
341 [arXiv: 100 1.2561 [hcp-ph]]. 

[23] B. Grzadkowski, M. Lindner and S. Theisen, Phys. Lett. B 198 (1987) 
64. 

[24] P. M. Ferreira and J. P. Silva, arXiv:1012.2874 [hep-ph]. 



22 



[25] S. Weinberg, Phys. Rev. Lett. 43 (1979) 1566. 

[26] W. Grimus and L. Lavoura, Eur. Phys. J. C 39 (2005) 219 [arXiv:hep- 
ph/0409231]. 

[27] P. Minkowski, Phys. Lett. B 67 (1977) 421. 

[28] T. Yanagida, in Proc. of the Workshop on Unified Theory and Baryon 
Number in the Universe, KEK, March 1979. 

[29] S. L. Glashow, in "Quarks and Leptons", Cargese, ed. M. Levy et al.. 
Plenum, 1980 New York, p. 707. 

[30] M. Gell-Mann, P. Ramond and R. Slansky, in Supergravity, Stony 
Brook, Sept 1979. 

[31] R. N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44 (1980) 912. 

[32] G. C. Branco, T. Morozumi, B. M. Nobre and M. N. Rebelo, NucL 
Phys. B 61 7 (2001) 475 [arXiv:hep-ph/0107164]. 

[33] J. A. Casas and A. Ibarra, Nucl. Phys. B 618 (2001) 171 [arXiv:hep- 
ph/0103065]. 

[34] M. Fukugita and T. Yanagida, Phys. Lett. B 174 (1986) 45. 

[35] F. R. Khnkhamer and N. S. Manton, Phys. Rev. D 30 (1984) 2212. 

[36] V. A. Kuzmin, V. A. Rubakov and M. E. Shaposhnikov, Phys. Lett. B 
155 (1985) 36. 

[37] J. Bernabeu, G. C. Branco and M. Gronau, Phys. Lett. B 169 (1986) 
243. 

[38] G. C. Branco, M. N. Rebelo and J. I. Silva-Marcos, Phys. Lett. B 614 
(2005) 187 [arXiv:hep-ph/0502118]. 

[39] F. del Aguila and J. Cortes, Phys. Lett. B 156 (1985) 243; G. C. Branco 
and L. Lavoura, Nucl. Phys. B 278 (1986) 738; F. del Aguila, 
M. K. Chase and J. Cortes, Nucl. Phys. B 271 (1986) 61; Y. Nir and 
D. J. Silverman, Phys. Rev. D 42 (1990) 1477; D. Silverman, Phys. 
Rev. D 45 (1992) 1800; G. C. Branco, T. Morozumi, P. A. Parada 



23 



and M. N. Rebelo, Phys. Rev. D 48 (1993) 1167; V. D. Barger, 
M. S. Berger and R. J. N. Phillips, Phys. Rev. D 52 (1995) 1663 
[arXiv:hep-ph/9503204]; M. Gronau and D. London, Phys. Rev. D 
55 (1997) 2845 [arXiv:hep-ph/9608430]; F. del Aguila, J. A. Aguilar- 
Saavedra and G. C. Branco, Nucl. Phys. B 510 (1998) 39 [arXivrhep- 
ph/9703410]; G. Barenboim, F. J. Botella, G. C. Branco and O. Vives, 
Phys. Lett. B 422 (1998) 277 [arXiv:hep-ph/9709369]; G. Barenboim, 
F. J. Botella and O. Vives, Phys. Rev. D 64 (2001) 015007 [arXiv:hep- 
ph/0012197]; G. Barenboim, F. J. Botella and O. Vives, Nucl. Phys. B 
613 (2001) 285 [arXiv:hep-ph/0105306]; K. Higuchi and K. Yamamoto, 
arXiv:0911.1175 [hep-ph]. 

[40] To appear in a future publication 



24 



